Parametrization of Nambu vertex in a singlet superconductor 

Andreas Eberlein^ and Walter Metzner^ 
^ Max- Planck- Institute for Solid State Research, D-70569 Stuttgart, Germany 

We analyze general properties of the effective Nambu two-particle vertex and its renor- 
malization group flow in a spin-singlet superconductor. In a fully spin-rotation invariant 
form the Nambu vertex can be expressed by only three distinct components. Solving exactly 
the flow of a mean-field model with reduced BCS and forward scattering interactions, we 
gain insight into the singularities in the momentum and energy dependences of the vertex 
at and below the critical energy scale for superconductivity. Using a decomposition of the 
vertex in various interaction channels, we manage to isolate singular momentum and en- 
ergy dependences in only one momentum and energy variable for each term, such that the 
OA singularities can be efficiently parametrized. 
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§1. Introduction 



In the last decade the functional renornialization group (fRG) has been estab- 
lished as a valuable source of new approximation schemes for interacting Fermi sys- 
o ■ tems.'^ These approximations are obtained as truncations of an exact functional 

^ . flow equation which yields the flow from the bare microscopic action to the quantum 

effective action as a function of a decreasing energy scale A^ 
C^ ■ Most interacting Fermi systems undergo a spontaneous symmetry breaking at 

sufficiently low temperatures (sometimes only at T = 0) . In the functional flow equa- 
tion, the common types of spontaneous symmetry breaking such as superconductivity 
^ ■ or magnetic order are associated with a divergence of the effective two-particle in- 

O , teraction at a finite scale Ac, in a certain momentum channel.'^^^ To continue the 

flow below the scale Ac, an appropriate order parameter has to be introduced. 

One possibility is to decouple the interaction by a bosonic order parameter field, 
via a Hubbard-Stratonovich transformation, and to study the coupled flow of the 
fermionic and order parameter fields. In this way order parameter fluctuations and 
^N ■ also their interactions can be conveniently treated. This route to symmetry breaking 

in the fRG framework has been explored already in several works on antiferromag- 
^ netisnP' and super conductivity.'^"'^ Typically the bare microscopic interaction can 

^^ ■ be decoupled in various channels, but the results obtained from truncated flow equa- 

^^ , tions depend on the choice of the Hubbard-Stratonovich field. This ambiguity is 

particularly serious in the case of competing instabilities with distinct order param- 
eters. 

k> , It is therefore worthwhile to explore also purely fermionic flows in the symmetry- 

5h I broken phase. This can be done by adding an infinitesimal symmetry breaking term 

to the bare action, which is promoted to a finite order parameter below the scale 
ylc.'^ A relatively simple one- loop truncation of the exact fRG flow equation was 
shown to yield an exact description of symmetry breaking for mean-field models such 
as the reduced BCS model, although the effective two-particle interactions diverge at 
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the critical scale ylc.'^''^ It turned out that the same truncation, with a very simple 
parametrization of the effective two-particle vertex, provides also surprisingly accu- 
rate results for the superconducting gap in the ground state of the two-dimensional 
attractive Hubbard model at weak coupling.!^ 

The purpose of the present work is to further develop the fermionic fRG route 
to spontaneous symmetry breaking, focusing on the case of singlet superconductiv- 
ity as a prototype for a broken continuous symmetry. We stay with the one-loop 
truncation used already in the previous works, but we improve on the parametriza- 
tion of the effective two-particle interaction in two directions. First, we make full 
use of spin-rotation invariance to reduce the number of independent Nambu com- 
ponents of the two-particle vertex to a minimum. Second, we pave the way for an 
efficient parametrization of singular momentum and energy dependences by extend- 
ing a recently proposed channel decomposition^^' of the two-particle vertex to the 
symmetry-broken state of a superconductor. We gain insight into the singularity 
structure of the vertex by analysing the exact fRG flow of an extended mean-field 
model featuring not only reduced BCS interactions, but also charge and spin forward 
scattering. 

The paper is structured as follows. In §2 we review the one-loop truncation of the 
exact flow equation, and we introduce the Nambu representation for the two-particle 
vertex in a spin-singlet superconductor. Available symmetries are exploited in §3, 
with the aim of parametrizing the Nambu two-particle vertex by a minimal number 
of independent components. In §4 we compute the exact renormalization group 
flow for the reduced model with BCS and forward scattering interactions, paying 
particular attention to the singularities of the Nambu vertex developed in the course 
of the flow. The channel decomposition introduced by Husemann and Salmhofer^ 
for an efficient parametrization of momentum and frequency dependences in the 
normal state is extended to singlet superconductors in §5. We finally conclude with 
a summary of our results in §6. 

§2. Truncated flow equations 

We consider a system of interacting spin-g fermions with a single-particle basis 
given by states with a momentum k, a spin orientation a G {tii}) ^^^ ^ kinetic 
energy e^. The system is specified by an action 

S[i^, i^] = ^{-iko + a) ^kai^ka + V[i^, i;] , (2-1) 



where k = (/cq, k) contains the Matsubara frequency ko in addition to the momentum; 
^|Jka and ipka are Grassmann variables associated with creation and annihilation 
operators, respectively; ^k = ^k — /^ is the single-particle energy relative to the 
chemical potential, and y[V', "0] is a two-particle interaction of the form 

^[V',V'] = ^ J2 J2 ^^I<^2'73a4{kl,k2,k3,k4)i^kiai'4'k2a2i^k3a3i^k4a4 ■ (2-2) 

ki,...,k4 o-i,...,CT4 
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Here and in the following the usual temperature and volume factors are incorporated 
in the summation symbols. 

The starting point for our analysis of the interacting Fermi system is an exact 
flow equatioiF^''^ for the effective action F [^/J, ip], that is, the generating functional 
for one-particle irreducible vertex functions in the presence of an infrared cutoff A. 
The cutoff is implemented by endowing the bare propagator with a regulator func- 
tion, r [Tpjip] interpolates smoothly between the bare action F °[ip,ip] = S[ip,ip] 
and the final effective action r[ip,ip] in the limit yl — )• 0. Spontaneous symmetry 
breaking can be treated by adding an infinitesimal symmetry breaking field to the 
bare action, which is then promoted to a finite order parameter in the course of the 
fiow.GO 

Expanding the exact functional fiow equation for F [ip, tp] in powers of ip and ip, 
one obtains a hierarchy of fiow equations for the n-particle vertex functions /^(2")^, 
We truncate the hierarchy at the two-particle level, including however self-energy 
corrections due to contractions of three-particle terms.l^ This truncation, which 
is sketched diagrammatically in Figs. 1 and 2, was used in all previous studies of 
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Fig. 1. Diagrammatic representation of the flow equation for the self-energy. The slashed line 
represents the single-scale propagator 5* . 
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Fig. 2. Diagrammatic representation of the flow equation for the two-particle vertex. The dots near 
the internal lines indicate differentiation of the propagator products with respect to the scale A. 

symmetry breaking from a fermionic flow.'^"'^ The fiow of the self-energy E is 
determined by a tadpole contraction of the two-particle vertex F^'^' and the socalled 
single-scale propagator 

S = OaG I17/1 fixed) (2-3) 

where G is the infrared regularized propagator. The latter is defined by (G )^^ = 
F^"^' -irR with a suitable regulator function R . The fiow of the two-particle vertex 
pW^ is determined by a one-loop diagram involving F^^' itself and the total cutoff 
derivative of the regularized propagator 
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The truncated system of flow equations described in Figs. 1 and 2 has the merit 
of solving mean-field models such as the reduced BCS modeP^J and the reduced 
charged density wave modeP^i exactly. 

We now specify to singlet superconductivity, where the continuous U{1) symme- 
try associated with charge conservation is spontaneously broken, while spin-rotation 
invariance remains conserved. It is useful to use Nambu spinors (f)ks and (pks instead 
of tl^ka and f/^fco- as a basis. The two basis sets are related by 

4>k+ = V'fct) 4>k+ = V'fet' 4>k- = 1p-ki, <t>k- = i^-kl ■ (2-5) 

The effective action as a functional of the Nambu fields, truncated beyond two- 
particle terms, has the form 



r^icp, ^] = r(o)^ + Y.Y. 4'i^(^) 4..<^fcs. 

k si,S2 

+ 4 X] X] ^stlts3sA^l^k2,h,k4,)(j}k^si4>k2S24>k3S34'k4Si , (2-6) 

fci,...,fc4 Si,...,S4 

where F^-^' yields the grandcanonical potential. Due to spin-rotation invariance only 
terms with an equal number of (p and (p fields contribute. The regularized Nambu 
propagator G^ is related to T^"^'^ by (G^)~^ = r*-^)"^ + R^, and can be written as 
a 2 X 2 matrix of the form 

C^(.^ _ ( Gi+ik) Gi_{k) \ _ ( G\k) F\k) \ 

^ ^^>-\ G^^{k) G'^_{k) ) ~\ F*^{k) -G^(-k) ) ■ ^^''> 

Note that the anomalous propagator F (k) is a symmetric function of k^ and k, due 
to spin-rotation and reflection invariance. 

For a singlet superconductor, the flow equation for the Nambu self-energy SI, 
shown graphically in Fig. 1, is given by 

"7T^sis2V^) — ~ 2^ 2^ ^s4s-3\'^ )-'siS3S4S2(^' k ,k ,k) , (2-8) 

k S3,S4 

and the flow equation for the Nambu vertex (Fig. 2) reads 

"^ rJtltsssS'^u k2, ks, k,) = il™, ,(A:i, A;2, k^, k,) - <f,',^r,,(A:i, ^2, fcg, h) 



dA 

~ 2 



1 _pp 



7;n^,s2S3s,{kiMMM) , (2-9) 



where 



<f.t4(^i>^2,fc3,A;4)=E E J-A^G\s'SP)Gi4^)] 



dA' 

X ^i'44s4(^l'P'9,fc4)r!Ji,,(^,A:2,fc3,p) , (2-10) 
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(2-11) 



<f..3.4(^i>^2,fc3,A;4) = E E ;^[^f;4(^)Gi4(«)] 



P'9 4, ■■■,4 



(2-121 



Note that the particle-particle terms in the Nambu representation contain particle- 
hole contributions in the original fermion basis and vice versa. In particular, the 
original particle-particle contribution driving the superconducting instability is con- 
tained in the Nambu particle- hole diagrams. 

We assume translation invariance such that the one-particle propagator and 
self-energy depend only on a single energy and momentum variable k, while the 
two-particle quantities are non-zero only if A;i + A:2 = ^3 + ^4, and can therefore be 
parametrized by three independent energy and momentum variables. 

§3. Constraints from symmetries 

In this section we exploit the available symmetries to reduce the number of 
independent components of the two-particle vertex in the superconducting phase to 
a minimum. In any case, the effective action is invariant under translations, spin 
rotations and spatial inversions. In most singlet superconductors, also time reversal 
invariance remains unbroken. 

In addition to the normal interaction involving two creation and two annihila- 
tion operators (^■0'0V) there are also anomalous vertices corresponding to operator 
products ijjTp-ip'tp + conjugate and -tpipTpip + conjugate.'^''^ We now write down 
manifestly spin-rotation invariant forms for each of these terms in the ^/'-basis, and 
then translate to the Nambu representation. 

A spin-rotation invariant normal interaction can be written as^ 

r(2+2)[V',^] = 7 E i^i^i^ ^2, k3, k4,)5a^aja2a3 " ^(A;i , /C2, ^4, A;3)(5<^ia3<Ja2<74] 

X ^k^a3k2U2'^hia-i'>Pkiai ■ (3-1) 

Here and in the remainder of this section we suppress the superscript A denoting 
the scale dependence. Alternatively one may write r^'^~^^[ilj,'ip] as a sum of a spin 
singlet and a spin triplet components' 



r(^+^)['0,V'] = 7 E [^'^(^i'^2,fe3,A;4)5o-i^2o-3c.4 + ^^(A:i,/c2,A:3,/c4)r^i ^2^.30-4] 

X V'fei(7iV'fc2a2'0A:3<T3V'fc4(T4 ' (3-2) 
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and 

V'^{ki,k2, h, h) = V{ki, k2, ks, k4) + V{ki, k2, k4, ks) , (3-3) 

V^{ki,k2, ks, ki) = V{ki, k2, ks, k4) - V{ki, k2, k4, ks) . (3-4) 

Time reversal invariance^ and charge conjugation (corresponding to hermiticity of 
the underlying Hamiltionian) yield the following relations for V{ki,k2,k^,k4): 

V{ki,k2,k3,k4) = V{kJ ,kj ,k2 ,ki) (time reversal) , (3'5) 

V{ki,k2,k3,k4) = V*{kl,k^,k2,kl) (charge conjugation) , (3-6) 

where k^ = (/cq,— k) and k* = (— /co,k), and the same for V^ and V'^ . Note 
that spatial inversions also transform k to k', but without permuting the momenta 
in V{ki,k2,k3,k4). The invariance of the action under the exchange of identical 
particles yields separate symmetry relations for V^ and V'^ under ki o k2 and 

^3 •(-)• k^: 

V^ikuk2, ks, k4) = V^{k2, ki,k3, ki) = V^{ki,k2, ki, ks) , (3-7) 

V^{ki,k2, ks, ki) = - V^{k2, ki,k3,k4) = -V^{ki,k2, ki, ks) , (3-8) 

while V obeys only 

V{ki,k2, k3, ki) = V{k2, ki,ki, ks) . (3-9) 

A spin-rotation invariant ansatz for the anomalous interactions can be con- 
structed systematically by requiring that commutators with the three components 
of the total spin operator vanish. 

A spin-rotation invariant anomalous interaction involving four creation or four 
annihilation operators can be written in the form 

^'^[W^{ki,k2,k3, A;4)(V'fcit^fe2; - i'kiii'k2t)i'4'ksti'k4l - V'fcgiV'feit) 

- w'^{ki,k2, ks, ki) [{'4>kit'>Pk2i + '4'kiii>k2t)(.i'k3^i>k4i + V'fc34V'fc4t) 

-2(V'fcitV'fc2tV'fc34.V'fc44 + ^fci;V'fc2i^fe3tV'fe4t)] +conj.| , (3-10) 

consisting of a singlet (S) and a triplet (T) term, which are separately invariant 
under spin rotations. Charge conjugated terms are denoted by "conj." . Time reversal 
invariance yields a relation between W ' and W ' , which assumes the simple form 

W^'^'^ku k2, k^, ki) = Ty^'^(-fci, -k2, -ks, -ki) , (3-11) 

if the order parameter (gap function) is chosen real. Furthermore, the functions 
W^{ki, k2,kz,ki) and W'^{ki, k2,k3,ki) are symmetric and antisymmetric under the 
exchange ki -H- k2 or A;3 o ki, respectively, and symmetric under the simultaneous 
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exchange ki -H- k^ and ^20/04: 

W'^iki, k2, k3, kA) = W^{k2, kiM, h) = w^ikiMM. k3) 

= W'^{k3,k4,ki,k2), (3-12) 

W^'iki, k2, ks, ki) = -W^{k2, ki, ks, ki) = -W^{kiM. h, ks) 

= W^{k3,k4,ki,k2). (3-13) 

Finally, a spin-rotation invariant anomalous interaction with three creation and 
one annihilation operators, or vice versa, can be written as 

r(^+^)[V',V'] = 2 X^l^'^^^l' ^2,^3,^:4) ^V'feia(V'fc2t^fc34-V'fc24V'fc3t)V'fc4^ 
ki a 



+ 2(V'fcltV'fc24,V'fc34V'fc44, - V'fcuV'fc2tV'fc3tV'fc4t) 



+ conj 



].}, (3-14) 



where e^ = 1 and e^ = —1. The terms with the coefficients X and X are separately 
spin-rotation invariant. Time reversal invariance yields the relation 

X^'^*(A:i, ^2, ^3,^4) = X^'^{-ki, -k2, -k3, -ki) (3-15) 

for a real choice of the order parameter. Invariance under particle exchange yields 



X'^^iki, k2, ks, ki) = X'^^iki, /C3, k2, ki) , 



X 



T t 



Ti 



,k2,k3,ki) = -X {ki,k3,k2,ki) . 



(3-16) 



To express the two-particle interaction in terms of Nambu fields, it is convenient 
to collect the 16 components of the Nambu vertex rsis2S3S4 in a 4 x 4 matrix 
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J- J L^ 
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(3-17) 



V r^4^_ r^4__ r['l^_ r^L_ J 



Note that the rows in this matrix are labelled by si and S4, while columns are 
labelled by S2 and S3. With this assignment the Bethe-Salpeter equation yielding 
the exact Nambu vertex in reduced (mean-field) models can be written as a matrix 
equation. Translating the spin-rotation invariant structure of the various interaction 
terms described above to the Nambu representation, one obtains a Nambu vertex of 
the following form 



T^^KkiMMM 
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( V^{h,k2,h,h) X{ki,k2,h,h) X*{klklklkl) -V{ki,h,h,k4) \ 
-X{ki,k2,k4„k3) W{ki,k2,k3,k4) V{ki,k4,k2,k3) X*{k4,k3,ki,k2) 
-X*{klkl,k*2,kt) V*{ki,h,k2,ks) W*{klklk*2,kt) X{kl,kl,kl,kl) 

\-V*{ki,h,k2,ki) -X*{ki,k3,k2,ki) -X{kl,klklkl) V^*{ki,k2,k3,ki)/ 

(3-18) 

where k = —k. The matrix elements W and X are related to the anomalous (4+0) 
and (3+1) interactions, respectively: 

W{ki,k2, hs, ki) = W^{ki, -ki, -k3, k2) - W^{ku -ks, -ki, k2) 
+ W^{ki, -k4, -ks, k2) - VF^(fcl, -fes, -^4, k2) 

+ 2W^{ki,k2,-k3,-ki), (3-19) 

X{ki,k2, ks, ki) = X^{ki,k2, -fca, k^) - X^{k2, fci, -/ca, k^) 
+ X'^{ki,k2, -ks, ki) - X'^{k2, ki, -ks, k^) 
+ 2X^{-k3,k2,ki,k4) . (3-20) 

The functions W and X obey the following relations under exchange of variables: 

W{ki,k2,k3,k4) = -W{k2,ki,k3,k4) = -W{ki,k2,k4,k3) 

= W{-k4,-k3,-k2,-ki), (3-21) 

X{ki,k2, k3, ki) = -X{k2, ki,k3, ki) . (3-22) 

In summary, by fully exploiting the available symmetries, in particular spin- 
rotation invariance, the number of independent functions parametrizing the Nambu 
vertex could be reduced substantially compared to what one gets by exploiting only 
the conservation of the z-component of the spin.'i^ 

§4. Reduced BCS plus forward scattering model 

In this section we consider a model with reduced interactions in the Cooper and 
forward scattering channels. The truncated flow equations described in §2 yield the 
exact flow for that model. By analyzing the flow we gain insight into the role of the 
various anomalous interactions and their singularities at and below the critical scale 
for superconductivity. 

The model is defined by the following action 

5[V',V^] = So[il^,^] + Vsc[^,^] + Vc[tP,i;] + Vs[tP,ij] , (4-1) 

consisting of a quadratic and three interaction terms. 

So[i^,i^] = ^(-iA:o + ek)4<x^fca + J^ (^(°Hfc)4tV^-H + ^^°^*WV'-H^fet) (4-2) 

k,cr k 
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contains the kinetic energy and an external pairing field A^^' , while Vsc, Vc and Vg are 
reduced interactions in the Cooper (superconducting), forward charge, and forward 
spin (magnetic) channel, respectively: 



Vsc[i^,i^] = ^ V"(°)(A;i,A;2)^fcitV'-fci4,V'-fc24.^fc2t ' (4-3) 

fci,fc2 

Vc[lp,1p] = ^J2 J2 -^i°Hfcl)^2)V'fciai^fc2<T2V'fc2^2V'fcl(7l , (4-4) 

ki,k2 o"i,(T2 

Vs[i^,i^] = 2Yj 12 Yj ^i°H^l'^2) (V<x; -^(72^^) '4^kiai-lpk2a2i^k2a'^i^kia[ , (4-5) 
kl,k2 criyCr'-^ 'T2,o-2 

where r = {t^ ^t^ , r^) is the vector formed by the three Pauli matrices. The coupling 

functions V^^>{ki,k2)i Fc {ki,k2) and Fs {ki,k2) are arbitrary functions of ki and 
k2- The interaction terms correspond to a spin-rotation invariant normal vertex of 
the form (13-11) with 



V{ki,k2,ks,k4,) = V^^\ki,kA)6k^^^k2h3-ki + F^^\ki,k2)5k^M^k2M 

- FP{ki,k2){Sk,,kJk2,k, + 26k,,k3^k2M) ■ (4-6) 

The pairing field A^^' breaks the U{1) charge symmetry explicitly. Spontaneous 
symmetry breaking is obtained in the limit A^^' — )• 0. The model (|4-ip is a gen- 
eralization of the reduced BCS model, where only Vsc contributes, whose flow was 
discussed extensively in Ref. 10). A special version of the model (j4-ip with Vg = 
was solved by a summation of all contributing Feynman diagrams in Ref. T8) . 



4.1. Exact integral equations and Ward identity 

Due to the restricted momentum dependence of the interaction terms, all contri- 
butions to F^"^' and F^'^' discarded in the truncation described in §2 vanish in the 
thermodynamic limit. This follows from a straightforward generalization of the argu- 
ments given in Ref. [T0|) . The constrained momentum dependence of the interactions 
leads to the following constraints on the external momenta of the various one-loop 
contributions to the flow of the Nambu vertex Fs^s2S3S4{ki, k2, k^, k/C), Eq. (|'2-9p : 

nsi^;%si{ki,k2,k-i,ki) oc Ski,k^Sk2,k3 , 
n^J^f^^^{ki,k2,k3,k4^) oc 5k^,k.j,5k2M > 
-^5^^25334(^1' ^2, k^, ki) oc 5k-,-k2hz-kA ■ (4-7) 

It is sufficient to consider the flow equation for the Nambu vertex with ki = k^ 
and k2 = k^, since the non- vanishing matrix elements for other choices of momenta 
follow from symmetries. Choosing ki = k^ =: k and /c2 = k^ =: k', only the direct 
particle-hole term contributes and the flow equation simplifies to 

'TJ^siS2S3S4\f^^ k ,k ,k) = ng_^g^g^g^{k, k ,k ,k) = 



10 



A. Eherlein and W. Metzner 



(4)A 



E J7[^44(^)^44(^)]^i.44^.^''^'^'^)4t.34(^'^''^''^) • 



(4-8) 



P'« 



This differential equation is equivalent to the Bethe-Salpeter-hke integral equation 

^sis2S3S4\'^^ k ,k ,k) = rg_^g^g^g^{k, k ,k ,k) 

summing up all Nambu particle- hole ladder diagrams. Note that G " = 0. Inserting 
this implicit solution for F^^' into the flow equation (|2-8|) for the self-energy, and 
using the relation S = G +G Jl G , the flow of the self-energy can be integrated, 
yielding 

Ei„^{k) - E^M = - E E Gis,ik')r^tlX,,ik, k', k', k) , (4-10) 

k' S3,S4 

with the bare Nambu vertex on the right hand side. This is just the familiar mean- 
field equation for the self-energy, which is exact for the reduced model ()4-ip . 

The exact solution determined by Eqs. (j4-9p and (j4T0p fulfils the Ward identity 
following from global charge conservation,^^' 

A^k) - zAW(fc) = ^Y1 [^^'Hk')G^+ik')GUk') - A^^>{k')Gt{k')GiAk') 

k' s,s' 

^^ +o'o_("') k , k , k) , 



(4-11) 



which connects the anomalous self-energy E^_{k) = A'^{k) with the two-particle 
vertex. The Ward identity implies that some components of the Nambu vertex 
diverge in case of spontaneous symmetry breaking (Z\ finite for A^^' — t- 0). 

4.2. Explicit solution and singularities of Nambu vertex 

An explicit solution of the Bethe-Salpeter equation (|4-9p for the Nambu vertex 
can be obtained for the case of separable momentum dependences of the interaction 
terms in the reduced model. In this special case the singularities of the various 
components of the Nambu vertex become particularly transparent. 

Using the 4x4 matrix representation of the Nambu vertex, Eq. (|3-17p . one can 
write the integral equation ()4-9p in matrix form as 

(4-12) 



r(4)^(/fc, k') = r(4)^« {k, k') + Y, r(')^°(A:,p) L^{p) r^^^^ip, k') 



where r(^)^(/i:, k') = r(^)^(A;, k' , k' , k) and L^ is the 4 x 4 matrix defined by 

C<A r^A r^A /-~<A r~iA r^A 

Lt Lj I 1 Lt I Lt I Lj Lt I 

r^A /-~fyl r^A r^A /~iA /— f/l /^A /~iA 

\ Lj I \J I Lt Lj 1 Lt I Lt Lj Lt / 



/ riA r'A 
r~<A r^A 
r^A /-^yl 



r^A /^A 
r<A r~<A 



r<A /-^yl 

r<A /-(/I 
r<A r<A 



(4-13) 



The bare Nainbu vertex r*^^-* " 
Eq. (|4-6p has the form 



{k,k') 
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k') corresponding to the reduced interaction in 



.{4)/lo 



where 



/ Fl' 


>{k, 



k') 













k') 


Ff 





k') \ 









y(0) 


*{k, 


k') 














[ Ff ) 


*{k 


k') 














^(0) 


*{k 


k') ) 



(4-14) 



F['\k,k') 
Fi'\k,k') 



FPik,k') + F(^Hk,k'), 
F(^)(k,-k')-F('\k,-k') 



(4-15) 



The structure of the full Nambu vertex T^'^' {k,k') is obtained by specializing the 
general form Eq. ()3-18p to the case ki = k4, = k and k2 = k^ = k' as 



r(4)^(A:,A:') 



/ F^'ik.k') 
X^{k',k) 
X^*{k'*,k*) 
V F2^*ik,k') 



X^{k,k') 

W'^{k,k') 

V^*{k,k') 

-X^*{k,k') 



X^*{k*,k'*) 

V'^{k,k') 
W'^*{k*,k'*) 
-X^{k*,k'*) 



Fi{k,k') \ 
-X^*{k',k) 
-X'^ik'^k* 

F^*{k,k') J 



(4-16) 



In the following we assume that the bare coupling functions V^^' , Fc , Fs , and 
the external pairing field A^^' are real and frequency independent. Then T^^' {k, k') 
and A'^^k) are also real and frequency independent, such that the complex conju- 
gation operations in Eq. (|4T6p can be omitted. The matrix products in the Bethe- 
Salpeter equation ()4-9p can then be simplified considerably by employing the follow- 
ing orthogonal transformation 



U 





( ' 








1 \ 


1 





1 


1 





/2 





-1 


1 







Ui 








1 / 



(4-17) 



The transformed vertex r('*)^ = Ur*^"^^ U has the form 
/ FHk,k') 



fW^(A:,A:') 



A'^{k,k') 











-^'^{k,k') 



\ 






-2X'^{k,k') 

with the "amplitude" (A) and "phase" (^) components 

A^(fe, k') = V^ik, k') + W^{k, k') 
<P'^{k, k') = V^{k, k') - W^ik, k') 





-2X^{k\k) 


F^{k,k') J 



(4-18) 



(4-19) 
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and the linear combinations in the forward scattering channel 

F^{k, k') = Ff{k, k') ± F^{k, k') . 
The transformed bare vertex r'^^ "(fe, k') has only diagonal entries: 



(4-20) 



.(0) 



f(4M«(A;,fcO 



{k,k') 








\ 





y(o)(fe,fe') 














-y(o)(fc,fc') 














F^'\k,k') ) 



(4-21) 



The frequency summed matrix L (p) = ^ L (p) also transforms to a simpler 



block matrix structure 



L^(P) 



/ lUp 












\ 



V 





Lf(p) L{^ip 

L^ip) 

L^(p) Llip) J 



(4-22) 



where 



Lf (p) = Y.iGi-{p)? + Gi^{p)G^_{p) 
L^(p) = J;[G^_(P)]' - Gi4p)G^__{p) 
L^{p) = -2j2Gi4p)GUp) , 

PO 

Ll{p) = Y.[Gi^{p)]' + [Gi_ip)]\ 

Pu 

Ll{p) = Y,[Gi^{p)f -[GUp)f . 



(4-23) 



Po 



For an explicit solution, we now assume that the momentum dependences of the 
bare interactions factorize: 

y(0)(fe,A:')=5^°)/sc(k)/.c(k'), 

FW(fe,fe')=#/e(k)/e(k'), 

FW(fc,fe')=#/s(k)/.(k'), (4-24) 

where /sc(k), /c(k), and /s(k) are arbitrary reflection invariant form factors. Then 
the momentum dependences of the flowing interactions also factorize in the form 
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F^{k,k') = 2g^Uk)Mk'), 

F^{k,k') = 2g^Mk)f,{k'), 

W^{k,k') = g:^fU^)fsc{k'), 

X^(A:,A:') = 5x/c(k)/.c(k'), 
and also the linear combinations V it W , 

^^(fc,A;') = <7f/.c(k)/.e(k'), 

^^(fc,A;')=4/.c(k)/.e(k'), 
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(4-25) 



(4-26) 



with ga = 9^+ aw and g^ = g"^ - g^. 

The momentum dependences on the right hand side of the integral equation 
()4-9p now factorize, and the momentum integration can be isolated in the following 
numbers: 

p 



''X 

It 
l^ 



p 

5]Lf_(p)/2(p), 
p 



(4-27) 



Inverting the resulting matrix equation, one obtains the explicit solution for the 
flowing couplings: 

,(0)^-l 



{9'c 



2li 



ya — 


2d^ 


0^- 


1 


y^ 


(ff(°))-^+^^ 


9^ = 


It 

2d^ 


,7^ 


(9^'^)-' -It 


9c — 


2d^ 


n^- 


1 



where 



{9^'^)-'-lt\\{29f^)-'-lt\-{lt)" 



(4-28) 
(4-29) 
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Note that g^ and g^ are coupled to other interactions only indirectly via the prop- 
agators, while g^, g^, and g^ are coupled also directly. 

The mean-field equation for the self-energy (|4-10p implies that the interaction 
driven part of the gap function, A (k) — A^^' (k), adopts the momentum dependence 
of the form factor for V^^\k, k'). Assuming also A^^'{k) = 6^^' fscO^)^ we can write 

A^'ik) = 6^fsc{k) . (4-30) 

We now discuss the behavior of the various components of the Nambu vertex 
as a function of the flow parameter A, especially the singularities near the critical 
scale Ac for spontaneous C/(l)-symmetry breaking, at which the vertex diverges for 
^W _^ Q y/q consider the usual case where superconductivity is the only instability 
of the system, that is, no instabilities driven by forward scattering occur (g^ and g^ 
remain finite for all A). 

For A > Ac and A^^' = the anomalous propagator F{k) vanishes, such that 
l^ = — Z4, l-^ = 0, and /^ = if. The anomalous components of the Nambu vertex 
W and X also vanish, and 

The critical scale Ac is the scale where 

lA. = _ ^ G^^{p)G^^{-p)fl{v) = {g^'Y' , (4-32) 

p 

such that V{k, k') diverges. 

For A < Ac and/or A^^> ^ 0, anomalous components appear. Using the relation 

F{k) = A{k)[G(k)G{-k) + F{k)F*{-k)] , (4-33) 

one can write the gap equation contained in Eq. ()4T0p in the form 

5A _ 5(0) ^ _^(0)5^/^ . (4.34) 

Inserting this into the solution for <P {k, k'), one finds 

^Hk,k') = -^V^^\k,k'). (4-35) 

Note that the forward scattering interactions and the anomalous (3-|-l)-components 
of the vertex do not affect this result. For A'^^' — )■ 0, one finds {k,k') — )■ 00 for 
any A < Ac- This is the divergence required by the Ward identity ()4Tip following 
from global charge conservation, and is associated with a massless Goldstone boson 
in the symmetry-broken state. 

In contrast to the phase component ^ of the vertex, the amplitude component 
A is regularized by the gap below the critical scale. Slightly below Ac, it behaves 
as 

A'^ik, k') oc — ^- . (4-36) 

§^ + 0[{5AY] 
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For A^^' — )■ 0, this component is therefore of order (5 ) ^ shghtly below Ac. The 
anomalous (3+l)-interaction behaves as 



X^{k,k') 



oc 



^ + o[{6^y] 



(4-37) 



For A^^^ — )• 0, it is therefore of order {5^)~^ slightly below Ac. 

We finally illustrate the above results by plotting the flow of the self-energy 
and the vertex components for a specific choice of model parameters and at zero 
temperature. The kinetic energy enters only via the density of states, which we 
choose as a constant D{^) = 1 with |,^| < 1. All form factors in the interaction 
terms are chosen as unity, fscO^) = /c(k) = /s(k) = 1, and the bare couplings as 
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(0) 



-0.5, g. 



(0) 



0.2, and g. 



(0) 



0. The qualitative features of the flow do not 



depend on the size of the couplings. The spin coupling has been chosen zero because 
the spin channel is decoupled from the rest. The flow is computed for a sharp cutoff 
acting on the bare kinetic energy, such that |^| > A. Since all bare interactions are 
momentum independent, also the flowing self-energy and the various components of 
the flowing Nambu vertex are momentum independent. 

In Fig. 3 we show the flow of the gap 6 for various values of the initial gap 6^^' . 
Note the sharp onset of 6 at the critical scale Ac for 6^^' = 0. This singularity is 
obviously smeared out for 6^^' > 0. The flow of the phase component of the Nambu 



0.3 



0.2 



0.1 



^A, 5(0) = 

6'^, 5(0) = 10-4 

^A 5(0) ^ 10-^ 



0.1 



0.2 
A 



0.3 



0.4 



Fig. 3. Flow of the superconducting gap 5 for various values of the bare gap i5(°'. 



vertex ^ is shown in the left panel of Fig. 4, and its final value at 71 = as a 
function of 6^^' in the right panel. Note that (p has the same shape as S , and that 
^A=o (;iiverges as 1/6^^^ for small (5^^', as described by Eq. (j4-35p . In Fig. 5 we plot 
the flow of the amplitude component of the Nambu vertex A and of the anomalous 
component related to (3-|-l)-interactions X . Both diverge at the critical scale Ac 
for 6^^' — >• 0, but decrease again for A < Ac- Note that X is much smaller than A 
for any A. 
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< 



-5000 



-10" 



-1.5 X 10" 



Si") = 10--' 
5(0) = lo-" 
SW = 10-5 



0.4 



< 

I 



10** 



10^' 



10" 



100 




10-** 0.001 



Fig. 4. Left: Flow of the phase component of the Nambu vertex <P for various values of S^'^' . Right: 
Phase component of the Nambu vertex at the end of the flow (A = 0) as a function of 5^^'' . Note 
that <P is momentum independent due to the choice of a constant form factor /sc(k) — 1. 





— a(") 


= 10 


.... gm 


= 10 


a(») 


= 10 



Fig. 5. Flow of the amplitude component A (left) of the Nambu vertex and of the anomalous 
component X corresponding to (3-|-l)-interactions (right) for various values of 5'"'. Note that 
A and X are momentum independent due to the choice /sc(k) = /c(k) = 1. 



§5. Channel decomposition 

When solving the flow equation for the two-particle vertex for a realistic (not 
mean field) model, one faces the problem that rather complex momentum and fre- 
quency dependences are generated by the flow. In a translation invariant system, 
the two-particle vertex depends on three independent momentum and frequency vari- 
ables. In previous works the frequency dependence was usually discarded, and the 
momentum dependence was crudely discretized.'2I'T3'Ii2J That procedure can be jus- 
tifled by power counting, as long as the flowing interactions remain weak and regular. 
However, in case of spontaneous symmetry breaking the two-particle vertex diverges 
in certain channels and singular momentum and frequency dependences develop. 

In this section we present a first step toward an efficient parametrization of the 
two-particle vertex, which captures the singular momentum and frequency depen- 
dences associated with the fiow into a superconducting phase. The singular momen- 
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turn dependence of the vertices originates from that of the loop integrals in the RG 
equations. Therefore, we decompose the vertices into several interaction channels, 
with the purpose of isolating the singular momentum dependence of each diagram. 
Distributing the diagrams among the channels according to their momentum depen- 
dence then leads to channel-decomposed flow equations. 

For the normal vertices, we follow the approach by Husemann and Salmhofer^ 
and add several general two-particle interaction terms for the different channels to 
the bare interaction. The latter is not distributed among the channels to avoid am- 
biguities. The interaction channels then describe corrections to the bare interaction 
from integrating out modes in the functional integral. We make the ansatz 



r(^+^U[^,^] = v[^,ij] + lY.5k, 



9 / V -i^i-\-k2,k3^k4 



^aia4^a2crs 



X V'fciai^fc2(72V'fc3f73^fc4Cr4 , (5-1) 

where ^[■i/',^'] is the bare interaction and the superscripts sc, c, and s refer to "su- 
perconducting", "charge", and "spin" channels, respectively. The momentum (and 
frequency) argument in parentheses is the total momentum of the interacting elec- 
trons for the superconducting channel and the momentum transfer for the charge and 
spin channels. These are the variables for which a singular dependence is expected. 
If the other momentum dependences remain regular, the above ansatz provides a 
good starting point for a parametrization of the various interaction channels as in- 
teractions mediated by a boson exchange. 

After symmetrization and translation to Nambu notation, one of the Nambu 
components representing the (2+2)-interaction reads 

ri^l't~if'i,k2, h, h) = v^'Hh, -h, -k2, ks) + ^"£(^3 - k2) + v^;\iki - k^) 

- '^V^lUS^i + k2) - VulUS^i - h) , (5-2) 

where k^ = ki + k2 — k^. Note that the total momentum of electrons in the super- 
conducting channel has transformed to a momentum transfer in the Nambu repre- 
sentation, while the momentum sum ki + k2 now appears in one of the spin-channel 
terms. 

The propagators in Eq. (|2T0p transport the following momenta through the 
diagrams, which they depend on singularly: 



^Sr^4(^i'^2,A:3,A:4) 

n^^Zs3S4ih,k2,k3,kA) 



momentum transfer k^ — /c2 
momentum transfer ^3 — fci 
total momentum ki + k2 ■ 



The three terms contributing to the flow of the two-particle vertex in Eq. (j2-9p are 
distributed among the various interaction channels according to their momentum 
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dependence. Hence we assign the direct Nambu particle-hole diagram to the super- 
conducting channel, such that 

dA^kuhM = n^^!f_{ki,k3 - q, ks, ki-q) , (5-3) 

where the dependence onq = k3—k2 is expected to be singular in the superconducting 
phase. The Nambu particle-particle diagram is assigned to the spin channel, 

^Kl(^) = \^^^+^{k^,p-k^,ks+p,-ks) , (5-4) 

with p = ki + k2- The crossed Nambu particle- hole diagram determines the flow of 
the remaining two terms in Eq. ()5-2p . that is, 



^^.:l(9') - ^VkiM) = -<-^-(^i> -k2 - l', h - q', -h) , (5-5) 

with q' = ks — ki. Using the above flow equation for V^' , this yields 

J^^kiM*^^'^ = ■^n^-+-{h,q-h,k2 + q',-k2)-n^^f_{ki,-k2-q',ki-q',-k2) . 

(5-6) 
We have thus obtained flow equations for all three contributions V^'^' , V^' , and 
V^' to the normal interaction. These channel decomposed flow equations for the 
normal interaction are equivalent to those derived by Husemann and Salmhofer^ 
in the ^-representation. 

For the channel decomposition of the anomalous vertices, we profit from the 
symmetries described in § 3. For the functions X ' {ki,k2,k2,,ki) parametrizing 
7^(3+1) j^^^j jj^ gq^ ()3T4p . we expect a singular dependence on the variable A;2 + ^3, 
which is the total momentum of the Cooper pair contained in i~'''^"'"^^['0,^]. We 
therefore write 

X^'^{ki, k2MM) = Xl'X{k2 + ks) , 

X^'^(A:i, k2, ks, k^) = Xlll^{k2 + ks) , (5-7) 

where k^ = ki + k2 + k^. As a Nambu component representing the (3+l)-interaction 
we choose 

r||'_+ (/ci, A;2, h, ki) = X^{ki,k2, k^, k^) 

= X',;%ik2 - ks) + Xl;%{k2 - ks) 

- ^Sl(^i - ^3) - <■! (fei - ks) + 2XT;,t,,(/ci + k2) , (5-8) 

where now ki, = ki + k2 — k^. Comparing the momentum dependences with those of 
the bubbles, we assign 



d 



'SA ., , ^ , yTA 



rPH,d 



dA [K,M{k2-k^) + Xil%{k2-k^)\=n^X'\{kiMMM + k2-k^) , (5-9) 
d 



dA L 



^S,A n. u\ , vT,A n, u\ _ ttPH^cv ^ 



K,k (^1 - ^3) + ^4':^, {ki - k3)\ = 77;^i^;(A;i, A;2, k^, ki + k2 - ks) , (5-10) 
^^Tfctfc,(A:i + A:2) = -^i7^^_+(A;i, A;2, ^3, fci + A;2 - A;3) . (5-11) 
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Eqs. ()5-9p and ()5-10p are equivalent, since 

Ks2sss, {ki MMM) = nfJ;Zs, (^2 ,h,h,h). (5-12) 

Solving for X ' and X ' , and relabeling variables, one obtains 

^^fi'.yfca^^) = ^++-'+(^1' ^2, h -q,ki+q) + -n^^_^{q - k2, /c2, -h,ki + q) , 

(5-13) 

^<'1(9) = -^<+-+(^ - k2, k2, -k,M + q) ■ (5-14) 

For the anomalous (4+0)-interactions we expect a singular dependence on the 
momentum of each Cooper pair contained in /^'^^'')[^,^], Eq. (|3-10p . which leads us 
to the ansatz 

W^^^ikiM. ks, ki) = ^fc^;,lfc3(A:i + k2) , 

W'^'^{ki,k2, k3, ki) = VFjl,3(fci + k2) , (5-15) 

where k4 = —ki — k2 — k^. The corresponding Nambu component becomes 

ri^l:L{ki, k2, ks, k^) = w^{ki,k2,k3, ki) 

= Ki(ks-k2) + W^;l{ks-k2) 

- <'l(^i - ^3) - <l(fci - ks) + 2W^;lik, + k2) , (5-16) 

where k/j^ = ki + k2 — k^. Comparing once again the momentum dependences with 
those of the bubbles, we assign 



"^ Wk;%,ik3-k2) + w,^;^^^{k3-k2) = n'l1'^_{kiMMM+k2-k3) , (5-17) 



dA 
_d_ 

dA 



rS,A n 7 \ , inrTA n. i. \ ttPH4 



rS,A /, , X . yT,A n. i. \ nPH.cr 



KltS^i - ^3) + KCkS^i - ^3)J = KXf-{ki. k2, ks, ki + k2-k3) , (5-18) 

^^J,t(^l +^2) = -^n^^__{ki,k2,k3,k, + k2-ks) , {5-19) 



where Eqs. ()5T7p and (j5-18p are equivalent due to exchange symmetry. Solving for 
ly ' and W ' , and relabeling variables, one obtains 

dA^klM^^^ =n++L_{ki,k2 -q,k2,ki - q) + -n^^__{ki,q - ki,k2,q - k2) , 

(5-20) 

^<l('?) = -^<+-(fci'9-^i>fc2,(?-A;2). (5-21) 

So far no approximation has been made in rewriting the flow equations. Conse- 
quently, they capture the exact flow of the reduced (mean-field) model discussed in 
§ 4. There, however, redundancies in the Nambu vertex for the reduced model could 
be exploited in order to construct the solution entirely in the momentum channel 
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ki = k/i, /c2 = A;3, where the flow is determined exclusively by the direct Nambu 
particle-hole diagram. 

The channel decomposition of the vertex and the flow equations provides a very 
useful starting point for an efficient approximate parametrization of the momentum 
and energy dependences being generated in the various channels in the course of the 
flow for models with generic interactions, as shown for the normal (not symmetry- 
broken) state of an interacting Fermi system by Husemann and SalmhoferP^J 

§6. Conclusions 

We have addressed the problem of finding an efficient parametrization for the 
effective two-particle vertex in a spin-singlet superconductor, with the perspective to 
solve functional renormalization group flow equations for interacting Fermi systems 
with a superconducting ground state. We have constructed a manifestly spin-rotation 
invariant form of the vertex, which reduces the number of independent Nambu com- 
ponents to only three functions {V , W, and X). By studying the exact flow of 
the vertex for a reduced (mean-field) model exhibiting superconductivity and also 
forward scattering, we have identified the singularities of the vertex associated with 
the superconducting instability. We have then expressed the vertex as a sum of 
various interaction channels where potential singularities are isolated in only one 
momentum and frequency variable in each channel, and derived the corresponding 
channel-decomposed flow equations. 

Our work paves the way for a controlled solution of the rather complex flow 
equation governing the flow of an interacting Fermi system into a superconducting 
phase. Since singular dependences generated by fermion loops have been isolated 
in only one momentum and energy variable per channel, one can parametrize these 
singularities by a relatively simple ansatz with a tractable number of parameters. 
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